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The NASA Constellation program has relied heavily on Computational Fluid Dynamics 
simulations for generating aerodynamic databases and design loads. The Orion Project 
focuses on the Orion Crew Module and the Orion Launch Abort Vehicle. NASA TetrUSS 
codes (GridTool/VGRID/USM3D) have been applied in a supporting role to the Crew 
Exploration Vehicle Aerosciences Project for investigating various aerodynamic sensitivities 
and supplementing the aerodynamic database. This paper provides an overview of the 
contributions from the TetrUSS team to the Project Orion Crew Module and Launch Abort 
Vehicle aerodynamics, along with selected examples to highlight the challenges encountered 
along the way. A brief description of geometries and tasks will be discussed followed by a 
description of the flow solution process that produced production level computational 
solutions. Four tasks conducted by the USM3D team will be discussed to show how USM3D 
provided aerodynamic data for inclusion in the Orion aero-database, contributed data for 
the build-up of aerodynamic uncertainties for the aero-database, and provided insight into 
the flow features about the Crew Module and the Launch Abort Vehicle. 


Nomenclature 

A - Area, sq. in. 

Ai - rolling moment coefficient normalized by the rolling moment coefficient at a = 0° 

A n - yawing moment coefficient normalized by the yawing moment coefficient at a = 0° 

A Y - side force coefficient normalized by the side force coefficient at a = 0° 

C A - axial force coefficient 
C D - drag coefficient 
C L - lift coefficient 
C m - pitching moment coefficient 
C N - normal force coefficient 
C P - pressure coefficient 

CT - total coefficient of thrust for the Abort Motors 
D - diameter of the Crew Module, inch 
h - grid cell height, inch 

log(R/R 0 ) - convergence parameter using the L2-norm solution residual of all flow variables 
M, Moo - Mach number, freestream Mach number 
MW - molecular weight 
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Ncell - number of grid cells in a volume grid 

P, p - pressure, psf 

q - free-stream dynamic pressure, psf 

R - gas constant 

^universal - universal gas constant 

Re D - Reynolds number based on the diameter of the Crew Module 
S re f - reference area, JtD 2 /4, in. 2 
T - temperature, °R 

TR - total thrust ratio for the Attitude Control Motors 
u - velocity component in the x axis, fps 

Y - total velocity, fps 

v - velocity component in the y axis, fps 
w - velocity component in the z axis, fps 
X or x - streamwise coordinate, inch 

Y or y - spanwise coordinate, inch 
Z or z - vertical coordinate, inch 
a - angle of attack, deg. 

(3 - sideslip angle, deg. 

A - difference between two values 
(|) - roll angle, deg. 

Y - ratio of specific heat for a gas 
p - density 

oo - free-stream condition 

Subscripts: 
air - air gas 
e - exit plane 
jet -jet gas 
t - throat plane 
tot - total conditions 

Selected Abbreviations: 

ACM - attitude control motor 
ALAS - alternate launch abort system 
AM - abort motor 
AR - analysis request 

ARSM - nonlinear Algebraic Reynolds Stress Model 

BPC - Boost Protection Cover 

BC - boundary condition 

CAP - CEV Aerosciences Project 

CEV - Crew Exploration Vehicle 

CFD Computational Fluid Dynamics 

CM - crew module 

DES - Detached Eddy Simulation 

IML - inner mold line 

ID AT - Integrated Design Analysis Team 

LAS - Launch Abort System 

LAV - Launch Abort Vehicle (Crew Module with Launch Abort System) 

PD - pressure dilation 

SA - Spalart-Allmaras one-equation turbulence model 
SCC - Sarkar compressibility correction 
SST - Menter Shear Stress Transport model 
TetrUSS - Tetrahedral Unstructured Software System 
WTT - wind tunnel testing 

k-8 - Jones-Launder linear two-equation turbulence model 
k-oo - Wilcox 1988 linear two-equation turbulence model 

2 

American Institute of Aeronautics and Astronautics 



I. Introduction 

NASA’s Constellation program was initiated to address the needs of U.S. access to space in the post Space 
Shuttle era with the mission to the moon, Mars and beyond. Unlike the Space Shuttle, where both crew and cargo 
are launched simultaneously, the proposed architecture for the Constellation program included two separate launch 
vehicles: the Ares I crew launch vehicle and the Ares V heavy lift cargo vehicle. The Ares-I was to be the launch 
vehicle for the Orion Crew Exploration Vehicle (CEV) shown in Figure 1. The CEV consisted of 4 components 
which were the Launch Abort Tower (LAT), the Crew Module (CM), the Service Module, and the Spacecraft 
Adapter. The CM combined with the LAT comprises the Launch Abort Vehicle (LAV). While the Ares-I program 
has been canceled (including Ares I and Ares V), work is continuing on the Orion CEV. The intention is to mate 
them with a future launch vehicle. 

The NASA Constellation program has relied heavily on Computational Fluid Dynamics (CFD) simulations for 
generating aerodynamic databases and design loads. The aerodynamic database development of the Orion CM and 
LAV is managed at NASA Johnson Spaceflight Center under the CEV Aerosciences Project (CAP) where multi- 
center coordination is provided for wind tunnel testing (WTT) and computational support. 

The NASA Tetrahedral Unstructured Software System (TetrUSS) tool set 1 has served in a supporting role to 
CAP for investigating various aerodynamic sensitivities, such as WTT sting interference effects, WTT-to-Flight 
Reynolds number scaling with Abort Motors (AM), AM thrust offset, and Abort-Motor/ Attitude-Control-Motor 
(AM/ACM) jet interactions. Each of these tasks required the delivery of hundreds of flow solutions with stringent 
data requirements and very tight time schedules. In the initial stages, TetrUSS, being a research tool, was not 
equipped for such demands and required urgent developments of new capabilities. The transformation of 
TetrUSS/USM3D from a research-oriented tool to a production-level tool in the midst of a fast-paced database 
generation campaign is described in Ref. 2. This paper will provide an overview of the contributions from the 
TetrUSS team to the Project Orion CEV and LAV aerodynamics, along with selected examples to highlight the 
challenges encountered along the way. 

The paper is organized with a brief description of the CM and LAV geometries in Section II, followed by a brief 
overview of the TetrUSS system and capabilities in Section III. Section IV provides an overview of analysis 
requests that were assigned to the TetrUSS team, followed by Section V that describes the processes developed for 
producing production level computational solutions. Section VI describes the data used for comparison with four 
sample activities presented in Section VII. The Orion LAV and its aerodynamics have been designated sensitive but 
unclassified (SBU) due to export control restrictions under the International Traffic in Arms Regulations (ITAR). 
Therefore, to comply with SBU rules, some details have been removed from plots and figures in this paper. 

II. Geometry Description 

Several variations of the Orion CM and LAV configurations have been examined by USM3D. Figure 2 presents 
images of the Apollo CM and the Orion Integrated Design Analysis Team (ID AT) CM configurations. The Apollo 
CM had an apex cover whereas the Orion CM configurations did not. Figure 3 presents images of the 605-054 LAV 
and 606F LAV configurations. The Launch Abort Vehicle (LAV) is comprised of the CM inside the Boost 
Protection Cover (CM BPC) component of the LAT. The shoulder groove (see Fig. 3(c)) on the CM Base is not 
present in the early LAV configurations. Each nozzle has a jet boundary condition (BC) plane located far enough 
downstream of the nozzle throat that the inner nozzle wall is nearly conical from the jet BC plane to the jet exit 
plane (see Figs. 3(d-e)). 


III. Computational Tools 

TetrUSS was developed during 1990s to provide a rapid aerodynamic analysis and design capability to applied 
aerodynamicists. The system consists of component software for setting up geometric surface definitions 
(GridTool), generating tetrahedral grids (VGRID), computing Euler and Navier-Stokes flow solutions (USM3D), 
and extracting meaningful information from analysis of results (SimpleView). An overview of TetrUSS capabilities 
in early 2000 is presented in Ref. 1, with more recent enhancements in Ref. 2. A commercial tool for adapting 
tetrahedral grids to the flow solutions was also utilized to provide a key capability for the current work. The salient 
features of the tools are described below. 

A. Grid Generation 

VGRID 3,4 is a tetrahedral grid generator based on the Advancing Front Method (AFM) for generation of surface 
triangles and ‘inviscid’ field cells, and the Advancing Layers Method (ALM) for generation of thin-layered ‘viscous’ 
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cells. Both techniques are based on marching processes in which tetrahedral cells grow from an initial front 
(triangular surface mesh) until the volume around the geometry is filled. Unlike the conventional AFM, which 
introduces cells into the field in a totally unstructured manner, the ALM generates organized layers of thin 
tetrahedral cells, one layer at a time, while maintaining the flexibility of AFM. Once the advancing front process is 
completed in VGRID, an additional post-processing step is required using POSTGRID to close any open pockets 
and to improve grid quality. VGRID input files are generated by an interactive geometry manipulation program, 
GridTool 4 5 . This graphics tool can import surface definitions from Initial Graphics Exchange Specification (IGES) 
files containing Non-Uniform Rational B-Spline (NURBS) surfaces and curves, as well as PLOT3D point definition 
files. GridTool is used to manipulate the geometry and to define necessary geometric surface patches and grid- 
spacing (source) parameters. It uses OpenGL for 3D graphics, and is available for Mac and Linux systems. The 
graphical interface is based on the Fast Light Toolkit. 

B. Solution Adaptive Grids 

Solution adaptive gridding was accomplished by the CraftTech CRISP CFD® code 6 . CRISP CFD® is a stand- 
alone tool for mesh modification and quality improvement of three-dimensional mixed-element unstructured 
meshes. Meshes comprised of tetrahedral, prismatic, and hexahedral regions may be readily modified to generate 
more accurate flow solutions through local refinement and coarsening. The adaption is driven by either an estimate 
of the solution error, or gradients of flow features. For tetrahedral cells as used in the present work, local refinement 
is achieved by means of a constrained Delaunay refinement algorithm combined with a circumcenter point 
placement strategy. Any inconsistency between the circumradius of a tetrahedron and some desired point spacing 
triggers the point insertion procedure. This iterative cell refinement is repeated until the cell circumradii are 
consistent with the prescribed point spacing. Coarsening of the tetrahedral region is also permitted through an edge 
collapse procedure. In regions where the grid is distorted or where solution errors are negligible, edges may be 
selected for removal. All cells incident to the deleted edge are removed from the mesh, the adjacent cells are 
redefined, and the two nodes of the edge are collapsed to a single vertex. Capability exists to modify the viscous 
layers of tetrahedral grids and underlying surface mesh. However, initial experience led to only adapting in the 
inviscid region of the grid at this time. 

C. Flow Solver 

USM3D 7,8 is a parallelized tetrahedral cell-centered, finite volume compressible Reynolds- Averaged Navier- 
Stokes (RANS) flow solver. The term “cell centered” means that the finite volume flow solution is solved at the 
centroid of each tetrahedral cell. Inviscid flux quantities are computed across each tetrahedral cell face using various 
upwind schemes. Spatial discretization is accomplished by a novel reconstruction process, based on an analytical 
formulation for computing solution gradients within tetrahedral cells. The solution can be advanced in time by a 2 nd - 
order “physical” time step scheme, a 2 nd -order “dual” time step scheme, or to a steady-state condition by an implicit 
backward-Euler scheme. Several turbulence models are available: the one-equation Spalart-Allmaras (SA) model, 
and the linear two-equation k-e model by Jones and Launder, the Wilcox 1988 linear two-equation k-w model, the 
Menter Shear Stress Transport (SST) model, and the two nonlinear Algebraic Reynolds Stress Models (ARSM) of 
Girimaji and Shih/Zhu/Lumley. Detached Eddy Simulation (DES) has been implemented in all of the turbulence 
models. A capability to trip the flow at specified locations on aerodynamic surfaces has been implemented for the k- 
8 turbulence model, but fully turbulent flow is assumed for the results in this paper. USM3D also has capabilities for 
overset grids and dynamic grid motion, which are not utilized in the current work. 

IV. Summary of CAP Task Requests 

The requests by the CAP office at NASA Johnson Space Center for computational support from the TetrUSS 
team are chronicled in Table 1. The list contains an identifier for each task, the subject geometry, a descriptive title, 
and its main objective. The majority of the analysis requests (AR) was directed toward various configurations of the 
LAV. Table 2 contains the corresponding flow conditions, jet conditions, grid information, and date of the tasks. 
Note the majority of the data were generated at transonic and low supersonic speeds. The tasks covered multiple 
objectives such as adding data to the project’s aerodynamic database, calculating WTT sting interference 
increments, evaluating WTT-to-Flight Reynolds number scaling with the AMs firing, determining effects of the AM 
thrust offset, and investigating AM/ACM jet interactions. 
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V. Flow Solution Process 

The process for applying TetrUSS/USM3D toward generating accurate Navier-Stokes solutions with quantifiable 
uncertainty has evolved over the five years of supporting the Orion Project. Table 3 contains the flow solution 
process. The process is split into 3 stages: USM3D input stage (step 1), grid generation stage (steps 2-6), and 
convergence stage (steps 7-9). Steps 1-2 of Table 3 are conducted on a local computational platform whereas the 
remaining steps were conducted on the NASA Advanced Supercomputing (NAS) Columbia Altix cluster. Steps 4-8 
are controlled by shell scripts. The adaptive grid utility CRISP CFD® version 1.5.1 9 was run on a single processor 
while the remaining steps were conducted on 128 processors. Step 9 is completed manually. More details of the 
steps are provided below. 

A. USM3D Input Stage 

This section describes the current USM3D input strategy, which has evolved from several guideline studies and 
extensive application. Reference 10 contains more details on the various input files and techniques. 

1. Flow, Geometry, and Code Input 

Flow conditions (Mach number, Reynolds number, temperature, etc) and geometry conditions (i.e. angle of 
attack, reference area, etc) were specified in the ‘‘projectAnpC’ file for USM3D. Various code parameters, e.g. time 
step, flux scheme, under-relaxation factor, etc., were also controlled through this file. 

Flow re-initialization of selected user-prescribed volume cells was often used to alter default freestream 
conditions and enable more robust solution startup. This feature greatly improves the robustness of startup from 
freestream flow by re-initializing the cells within the high-pressure AM and ACM nozzles of the LAV, or in the 
subsonic pocket behind the CM Base embedded within a supersonic free-stream. 

During the guideline studies listed in Table 1, parametric studies were conducted on the available turbulence 
models in USM3D. For the tasks after Guideline Study 3 in Table 1, the SST turbulence model with the Sarkar 
compressibility correction (SCC) and pressure dilation (PD) was selected after statistical analysis of correlations 
with experimental wind tunnel data 10 . All computations were performed with the assumption of fully turbulent flow. 
Because of the very tight time schedules, virtually all runs were made “steady-state” with local time-stepping 
convergence acceleration. An occasional time-accurate solution was conducted when unsteady flow effects were 
significant. 

2. General Boundary Conditions 

The boundary conditions for the geometry surfaces and outer grid boundary were set in the (( project.mapbc ,, 
input file that lists all of the surfaces of the grid. The outer boundaries of the grid were set to characteristic 
inflow/outflow condition. The surface boundaries of the vehicle were set to a no-slip BC except for the inner-mold- 
line (IML) surfaces of the ACM and AM nozzles. For jets-off conditions, the jet BC plane and the IML nozzle walls 
were set to a viscous boundary condition. For jets-on conditions, the jet BC plane has a boundary condition 
described in the next section and the IML nozzle walls were set to an inviscid boundary condition. If present, the 
sting surface was set to an inviscid tangent flow boundary condition. 

3. Jet Inflow Boundary Conditions 

The jet BCs for the AM and ACM nozzles were derived using the wind tunnel flow conditions or flight LAV 
trajectory data depending on the configuration and objective. Mach number (M), pressure (P), and density (p) were 
computed at the nozzle exits, then extrapolated using isentropic relationships to a jet BC plane located slightly 
downstream of the nozzle throat using isentropic relationships. The Appendix describes the equations and process 
used to calculate the jet BCs. The jet BCs were input to USM3D through the “project.jetbc” file. The velocity 
vectors on the recessed jet were “fanned” across the surface from the center nozzle axis until becoming tangent at 
the nozzle walls. No attempt was made to simulate the boundary layer of the jet inflow, thus a tangent flow 
“inviscid” BC was applied to the nozzle IML to improve numerical stability. 

The current USM3D code is limited to ideal gas with air, without a capability for modeling multiple species 
gases to simulate the chemical composition of the AM and ACM plumes characteristic of the flight vehicle. Hence 
the hot jet simulations are constrained to air with the ratio of specific heat of y a i r =1.4 and a gas constant of R a i r = 
1716.47. This assumption is reasonable for comparing with cold-jet wind-tunnel nozzle configurations. For the hot- 
jet USM3D simulations, the jet BC values for the non-air gas conditions( i.e. y jet and R jet ) are corrected from a non- 
air gas to air as described in the Appendix. The advantage of this type of correction for y and R is that it duplicates 
the true jet gas exit pressure - an important parameter in jet exhaust expansion and therefore plume shape. It also 
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duplicates the jet mass flow and exit velocity, both important parameters for entrainment and for maintaining desired 
levels of the thrust. 

B. Grid Generation Stage 

The grid generation process continually evolved as configuration geometries became more complex and high- 
powered jet plumes added new requirements for grid resolution. The following describes the range of techniques 
employed to deliver the required solutions on schedule. 

7. Base Grids 

The base tetrahedral grids (i.e. unadapted grids) were generated using the VGRID4.0 code 4 . Several configura- 
tions examined included wind tunnel stings. However, wind tunnel walls were never modeled. Full configuration 
tetrahedral grids were produced with an outer boundary box prescribed at minimum x, y, z of, at least -100D, -100D, 
-100D, respectively, and maximum x, y, z of, at least, 100D, 100D, and 100D, respectively, where D is the CM 
diameter. The grids were verified to have an average y-plus of 0.25 at the centroid of the surface tetrahedra. This 
parameter was mostly dependent on Re D . Thus for a common geometry, wind tunnel Re D conditions would have a 
different grid than those of flight Re D conditions. In the first year of the USM3D work, the base grid was used for 
the unpowered LAV configurations. 

2. Predefined Wake Grids 

The initial computations on a base grid for the Apollo CM exposed the critical need for grid refinement of the 
grid in the wake region. 11 This observation presented a considerable operational challenge since a large number of 
angles of attack (a) were computed between 0° and 180° with each requiring a unique grid refinement. To solve this 
problem, an h-refmement utility was applied in a four-staged process to progressively refine smaller and smaller 
sub-regions as illustrated in Fig. 4. Figure 5 presents several grids for the Apollo CM of Guideline Study l 11 of 
Table 1. Figure 5(a) presents the four successive grids generated in the predefined wake grid generation process for 
a = 0°. The directional refinement observed in Figs. 5(b-c) was generated automatically using a script for a = 40° 
and 150°, respectively. The script performs the proper user-prescribed rotations and successive subdivisions of the 
grid within the cone regions. 

3. Solution Adapted Grids 

As the CAP tasks became more focused on the LAV, grids became more difficult to generate due to increased 
geometric complexity and the challenge of volumetric grid refinement for components such as the AM and ACM 
plumes, and the cross-flow wake behind the tower and CM BPC. Several of the earlier tasks for the unpowered LAV 
(up to AR-60 in Table 2) were conducted with no adaption of the base grid. The eventual inclusion of powered jets, 
beginning with AR-76 12 , drove the firm requirement for solution adaptive grids to resolve the features and 
interactions of the AM and ACM plumes as illustrated in Fig. 6 for the 606F LAV configuration. Figure 6 presents 
results for the preliminary solution on the base grid at M^ = 1.1, a = 5.0°, Re D = 64.4 x 10 6 , a combined AM thrust 
coefficient (CT) of 4.0 for all 4 AM nozzles, and a combined thrust ratio (TR) of 0.04 for the two south-firing ACM 
nozzles # . Figure 6 shows data on a streamwise plane cutting through the AM nozzle that is offset from the symmetry 
plane by a roll angle (()>) of 45°. The grids and Mach contour data (on the § = 45° plane) are shown around the 
shaded LAV surfaces. The base grid and the Mach contour data from the preliminary solution are shown in Fig. 
6(a). The upper value of Mach number is clipped to a less than maximum level to enhance the details of the plume 
interactions with the aerodynamic surfaces, which results in the inner portion of the jet plume being solid red. The 
adapted grid based on the preliminary solution is shown in Fig. 6(b). The adapted grid is shown on top of the Mach 
contour data on the right side of Fig. 6(b). The clustering of adapted cells is evident in the AM and ACM plumes, 
around the shoulder of the CM BPC, and within the wake region of the CM Base. 

More detail of the 606F LAV grid is shown in Fig. 7. A layer of viscous cells near the surface exhibits large 
disparities in cell size relative to the adjacent adapted cells. This is an artifact of the CRISP CFD® code that only 
creates isotropic tetrahedra for the adapted cells. CRISP CFD® has the capability of adapting through the viscous 
layers, but has shortcomings with the surface grid adaption. The effect of adapting through the boundary layer was 
investigated and was determined to be negligible to overall force and moment values. Note the near absence of 
viscous layers at the intersection of the lower tower and fairing. The thickness of this layer is based on the relative 
size of the surface triangle compared to the nearby inviscid spacing - the smaller the surface cell size the thinner the 
viscous layer. 


# A discussion of the ACM jet firing orientation can be found in Ref. 13. 
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The length (as defined by number of iterations) of the preliminary solution, which was used to adapt the base 
grid, has evolved. For AR-76 through AR-105 the preliminary solution was a l st -order base-grid solution of 
approximately 3,000 iterations. From AR-112 to AR-128 the preliminary solution was allowed to progress to a 2 nd - 
order base-grid solution of approximately 15000 iterations. Figure 8 shows the coefficient convergence history of C A 
for two solutions that varied by the length of the preliminary solution used to derive the adapted grid. The solutions 
are for the 606F LAV configuration at = 1.1, a = 0°, Re D 64.4 x 10 6 , a combined CT = 3.0 for the four AM 
nozzles, and a combined TR = 0.04 for the two south-firing ACM nozzles. C A has a discontinuity where the solution 
switches from 1 st order to 2 nd order and at the restart on the newly generated adapted grid. The solution on the left of 
Fig. 8 had a preliminary solution of 3000 iterations and was a l st -order solution (referred to hereafter as the 1=3000 
solution). The solution on the right of Fig. 8 had a preliminary solution of 15,000 iterations and was a second order 
solution (referred to hereafter as the 1=15000 solution). 

Figure 9 presents Mach contour and adapted grid data from the preliminary 1=3000 and 1=15000 solutions of 
Fig. 8. The Mach contour data of the 1=3000 solution are shown on the left side of Fig. 9 and the data of the 1=15000 
solution are shown on the right side of Fig. 9. Beneath each Mach contour is the grid that was adapted using the 
solution at that iteration. Three different cross-sectional views are shown in Fig. 9. The two solutions show 
differences in the Mach contour data such as: 

• AM plume sits on the shoulder of the CM BPC for the 1=3000 solution but has lifted off that shoulder for the 
1=15000 solution (see Fig. 9(a)) 

• a shock on the upper surface of the upper tower is further downstream for the 1=15000 solution than in the 
1=3000 solution (see Fig. 9(b)). 

• ACM plumes are more symmetric for the 1=3000 solution than for the 1=15000 solution (see Fig. 9(c)). 

The above differing flow structures yielded slightly different adapted grids as seen in the accompanying adapted 
grid images of Fig. 9. These slightly different grids produce different convergence behavior as seen in Fig. 8. The 
1=15000 solution has asymmetric ACM plumes but C A levels out better than that of the symmetric ACM plumes of 
the 1=3000 solution. 

C. Convergence Stage 

The aerodynamic force and moment coefficients were extracted as an average of the last 1000 solution iterations. 
Thus, quantitative characteristics of the evolving behavior for these coefficients as the flow solution develops 
provide a numerical approach to determining convergence. Once an adapted grid was generated and initialized with 
the base-grid flow solution, the convergence stage (steps 7-9 of Table 3) began. After an additional 8000 iterations 
on the adapted grid, an automated monitoring script checked the USM3D output in 15-minute intervals to determine 
if the convergence criteria listed in Table 4 were met. When two of the first three criteria were met, the solution was 
automatically stopped. The last step of the run process required the user to examine a solution summary chart, and 
other data if necessary. If the user decided that more iterations were necessary to further the convergence of the 
solution, the solution was restarted and step 8 in Table 3 was repeated for a user-defined number of iterations. In 
cases where the solution fails, deviation from the default process was required. More details about the trouble- 
shooting process are given in Ref. 10. 

Once the monitoring script of Step 8 in Table 3 stopped the solution, the user examined the solution using a 
solution summary chart. A representative solution summary chart is presented in Fig. 10 for the unpowered 
ALAS1 lr3b configuration at M* = 1.3, a = -8°, and Re D = 2.5 x 10 6 at solution iteration 9,398. The right side of the 
summary chart displays flow Mach contour data on a streamwise plane cutting through the AM nozzle that is offset 
from the symmetry plane by a roll angle of 45°. In the lower right comer of the chart is a plot of the surface pressure 
coefficient (C P ) contours with a longitudinal cutting plane through the AM nozzles depicting Mach contour data of 
the jets. The left side of the summary chart tracks the convergence history of the solution. The lower left plot in Fig. 
10 depicts the history of the order-of-magnitude drop of the L2 residual norm of the conservation equations (five 
flow variables) and the upper left plot depicts the histories of the axial force coefficient (C A ), normal force 
coefficient (C N ), and pitching moment coefficient (C m ) for the nose-to-base configuration. The IML and jet BC 
plane for each motor and the sting were excluded from integration of the three coefficients shown in Fig. 10. 

The solution shown in Fig. 10 meets convergence criteria 1 and 2 at iteration 9,398. However, the delta of the 
maximum and minimum of C A over the last 1000 iteration is greater than the tolerance of 0.004 and thus criterion 3 
is not met. A visual inspection of Fig. 10 shows that C A has not achieved convergence, although the percentage 
requirement of 1% variation in criteria 1 and 2 are met due to the relatively large magnitude C A . Therefore a 5 th 
criterion was imposed on such cases. The 5 th criterion required that the absolute value of the delta between 
coefficients averages at the two iteration intervals be less than 0.002. The solution of Fig. 10 was extended until this 
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5 th criterion was met at iteration 11,219 as shown in Fig. 11 where the Ca convergence history levels off to within 
tolerance. About half of the solutions met two of the first three criteria but still did not appear converged upon visual 
examination of the convergence history. These solutions were extended until the 5 th criterion was met. 

Five different classifications of convergence history were observed in the USM3D tasks. The five types 
presented in Fig. 12 are: normal, oscillation, small oscillation, irregular oscillation, and curve break convergence. 
The type of convergence could sometimes affect the final aerodynamic coefficient averages from a solution. For 
example, a solution with a regular periodic oscillation in one or more aerodynamic coefficients with respect to 
iteration was considered to have an oscillation convergence (see Fig. 12(b)). Usually, the oscillating coefficients had 
two or more cycles in the last 1 000 iterations such that the average taken over the last 1 000 iterations was deemed 
acceptable. However, if less than two cycles occurred during the last 1000 iterations, a new interval of iterations 
based on the last two cycles would be used to calculate the averages and to re-compute the convergence criteria. The 
irregular oscillation convergence had no regular cycle in the oscillations (see Fig. 12(d)). The small oscillation 
convergence type (see Fig. 12(c)) was similar to the oscillation convergence type with a much smaller oscillation 
magnitude. The irregular oscillation and small oscillation convergence types had a coefficient oscillation less than 
0.004 and thus retained the 1000 iteration interval for calculating averages and convergence criteria. The curve 
break convergence type in Fig. 12(e) occurred with a sudden shift in the coefficient history, usually when a major 
shift in the flow structure developed in the evolving solution. Just before the break in the coefficient vs. iteration 
curve the convergence appeared normal and very close to satisfying the convergence criteria. After a sudden shift in 
the coefficient history the solution was allowed to run further until the convergence criteria is met. The curve break 
convergence type is not a common type. For AR-105, 11 of the 230 cases performed had the curve break 
convergence and 10 of those were at =1.1 and 1.3. The 1 1 th case occurred at JVU = 2.0 

VI. Comparison Data 

The USM3D results are compared with Project Orion data obtained from three sources. The first is experimental 
data from the various wind tunnel tests within the project. These tests provided force and moment measurements as 
well as limited surface pressure data. In the latter stages of the project, some pressure sensitive paint data were 
obtained but are yet to be compared with USM3D results. The second source of comparison data is from the Orion 
Aerodynamic Databook 14 which contains data also called the Orion aero-database. The aero-database is the project’s 
best estimate for the aerodynamics vehicle shown in Fig. 1. The aero-database is based on available experimental 
and computational data and has evolved as additional data are developed (including the USM3D tasks). The third 
source of comparison data is computational results from the OVERFLOW 15 structured-grid code which is producing 
large amounts of data to the Orion project. 

VII. Sample Applications 

Over the course of the Orion CAP Project, TetrUSS/USM3D has provided an extensive amount of computational 
aerodynamic and loads data on two types of configurations - the solitary CM and the LAV (with and without the 
AMs and ACMs firing). Approximately 500 solutions were obtained for CM configurations, 400 solutions for un- 
powered LAV configurations, and 1,100 solutions for powered LAV configurations. The following describes some 
typical results from the TetrUSS/USM3D tasks and their impact to the project. 

A. Apollo CM Configuration 

An initial guideline study was conducted on the Apollo CM to develop the practices for generating the static 
aerodynamic coefficients on CM configurations at transonic and low supersonic M m (i.e., 0.7 < M^ < 2.4). The 
practices were determined through systematic parameter studies of grid resolution, turbulence models, flux schemes, 
limiters, etc. The turbulence models that were considered are the Spalart-Allmaras (SA) model, the linear k-e model, 
the Girimaji nonlinear k-e model, and the Menter Shear Stress Transport (SST) model. The effect of turbulence 
model is examined in Fig. 13 for the Apollo CM configuration at M* = 0.9, wind tunnel R D , and 0° < a < 180°. The 
aerodynamic coefficients were calculated as the average over the last 1000 iterations. As observed in Fig. 13, no 
clear guidance is apparent regarding the selection of a superior turbulence model for the Apollo CM aerodynamics, 
illustrating the complexity of the flow field around such a benign geometry. Models that correlate well at one end of 
the a-range do not correlate as well at the other end of the range. From this study, both the linear and non-linear k-e 
models were eliminated due to lack of robustness. Hence, the SA and SST models were retained for future 
calculations on the CM geometry. 


8 

American Institute of Aeronautics and Astronautics 



Fig. 14 presents the Mach contour data on the Y = 0” cutting plane for twelve angles of attack. These images 
provide insight into the dominant flow mechanisms that induce the aerodynamic behavior. For example, the lower 
and higher angles of attack are dominated by pure blunt-base flows, whereas the mid-range angles (60° < a < 80°) 
experience a strong influence from smooth round surface cross-body flow separation, which is very difficult to 
model. The sensitivities in the mid-range are even evident as large differences in the experimental data presented in 
Fig. 13. Thus, one should expect any numerical prediction of aerodynamics for 60° < a < 80° to have large 
uncertainties. 

The variations in coefficients due to the different turbulence models were used in developing aero-database 
uncertainties. The CM simulation guidelines developed here were used to produce a full set of solutions on the CEV 
CM that were incorporated into the Orion aero-database. These guidelines also provided a starting point for the LAV 
guidelines. 

B. ID AT CM Configuration 

The ID AT CM configuration shown in Fig. 2 was examined in AR-122 of Table 1 to generate data for the most 
recent CM aerodynamic database. A small study was performed initially to refine the application guidelines. This 
task required a combination of predefined wake grids (used in the above sample) with additional adaptation of the 
flow solution. A base grid was first generated with a predefined cell clustering in the wake region at a prescribed 
angle of attack, and then CRISP CFD® was applied to further adapt the grid to the flow solution. Each angle of 
attack had a uniquely defined wake region. Figure 15(a) presents the baseline base grid and the baseline adapted grid 
for the CM at M* = 0.5 and a = 150°. Preliminary solutions on the ID AT CM indicated a strong sensitivity to the 
grid resolution around the CM shoulders. Therefore, a finer base grid was generated with additional clustering 
around the shoulders (see Fig. 15(b)). 

The USM3D aerodynamic coefficients (C A , C N , and C m ) are correlated with OVERFLOW data and the ID AT 
CM aerodynamic database Version 0.55 in Fig. 16 for M* = 0.5 at 135° < a < 160°. OVERFLOW predicts a break 
in the C A vs. a curve between a = 150° and a = 155°. Both steady state and time-accurate USM3D solutions are 
presented in Fig. 16 for the baseline and finer grids. The steady-state USM3D solutions do not predict, with either 
grid, the break in the C A vs. a curve observed in the OVERFLOW time-accurate solution. Also the increment in C A 
between the baseline and finer grids is nearly constant over the ot-range. In contrast, the increment in C N and C m due 
to grid resolution is significant and increases with decreasing angle of attack for the a-range shown in Fig. 16. 
Outside of 135° < a < 160° range, USM3D data agrees better with the OVERFLOW data (not shown here). 

Figure 17 presents Mach contour and surface pressure coefficient (C P ) contour data at a = 150° for both the 
baseline and finer grids. The lower surface of the CM body has a lower C P for the baseline grid solution which 
produces a higher C N and lower C m at a = 150°. The data of Fig. 17 represent the last iteration of a slightly 
oscillatory convergence. An examination of Mach and C P contours at several prior iterations verified that the 
observed behavior is not significantly changed. The variation in the coefficients’ oscillating convergence histories is 
due to small increases and decreases in the low C P data observed above. 

The convergence histories for the baseline and finer grids are presented in Fig. 18 for a = 150°. The steady-state 
solution is evident prior to the notation of the time-accurate solution in the figure. The time-accurate solutions were 
started from an iteration where the steady-state solutions were considered converged according to the criteria in 
Table 4. C A for the time-accurate solutions oscillates by approximately ±0.02 whereas C N and C m settled to within 
±0.003. Figure 19(a) presents the Mach contour data on the Y = 0” plane for the fmer-grid time-accurate solution at 
four time steps representative of two cycles of C N oscillation. In contrast to the steady solution (shown in Fig. 17), 
the wake in the time-accurate solution is highly oscillatory. The corresponding impact on the surface C P contours on 
the conical shell of the CM is presented in Fig. 19(b) as viewed from the front of the CM. Significant variation of C P 
with time is observed and could easily alter the interplay of forces that result in the breaks in coefficients of the 
time-accurate solutions in Fig. 16. 

To summarize, USM3D predictions, at these conditions, are sensitive to grid resolution, are predicting an 
unsteady flow about the ID AT CM, and do not compare well with OVERFLOW or aero-database results in the 
range ofl35°<a< 160°. These results indicate a need for experimental data to better quantify the aerodynamics in 
range of 135° < a < 160°. The USM3D IDAT CM results were used in defining uncertainties for the aero-database. 

C. 19-AA LAV Configuration 

Task AR-41 16 was initiated to provide supporting computational aerodynamic analysis of perceived WTT data 
anomalies and sting interference effects for the 19-AA LAV wind tunnel model. The 19-AA LAV model was based 
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on the LAV 605-054 PA-1 geometry but had a slightly longer upper tower. The 19-AA LAV model was tested at the 
Boeing Polysonic facility in St. Louis, Missouri. In order to completely cover the desired a-range of 0° to 180°, five 
different sting configurations were tested as presented in Fig. 20. Each sting had a different a-range as noted in Fig. 
20 with overlaps occurring. No testing occurred at 130° < a < 140°. Grids were generated using VGRID with no 
grid adaption. A sixth configuration with no support sting was added to the USM3D task and is referred to as the 
free-air configuration. 

Aerodynamic data over the full range of angle-of-attack at = 0.5, constructed from measurements using a 
composite of the five sting configurations, are presented in Fig. 21. The experimental data for the 5 different stings 
were plotted as solid lines in different colors. For example, measurements at a = 30° were made with both the 15° 
sting (red) and the 50° sting (green). Similarly the data at a = 60° and 70° reflect dual measurements from the 50° 
sting (green) and the 77.5° sting (blue). It is observed for C A , C N , and C m that the trends of the sting effect on the 
coefficients are consistent in the sting overlap regions between the USM3D solutions and the wind tunnel data, 
although the magnitudes are not always comparable. Such trends for sting effect are not fully consistent with C m as 
evident at a = 60°, 70°, and 95°. Sting-interference corrections relative to the free-air configuration were provided to 
the CAP team and incorporated into the aero-database for unpowered LAV flight. 

D. 606F LAV Configuration with ACMs and AMs Firing 

USM3D has primarily served as a second computational data source to the CAP team for supplementing the 
primary data computed with the OVERFLOW code. The objective for AR-112 in Table 1 was to provide data for 
assessing jet interference effects from the AM and ACM plumes on the LAV Fairing and CM BPC. Midway through 
the task, an additional focus was added to examine a flow asymmetry issue that was discovered during a parallel 
OVERFLOW task. During the OVERFLOW study, abnormal asymmetries were discovered in the ACM plumes at 
certain flight conditions for cases with symmetric jet conditions and a near symmetric geometry. Figure 3(b) 
presents sketches of the 606F LAV configuration for which the LAV/Service Module umbilical cavity is offset to 
the left of the LAV centerline. 

A short grid study was conducted to determine the sensitivity of AM/ACM interactions to grid resolution. One 
result from this study was the use of a longer, second-order solution more fully converged for the grid adaption as 
was discussed previously for Fig. 9. Figure 22 presents another aspect of the grid study - the use of a predefined 
cone grid on the upper tower. The baseline base grid and the base grid with a predefined cone grid (referred to as the 
predefined-cone base grid) are presented at two cross-section cuts in Fig. 22. Note that the predefined cone grid (on 
the left of Fig. 22) on the upper tower stopped just ahead of the AMs. 

Figure 23 presents the adapted grids and Mach contour data generated for the baseline and predefmed-cone grids 
at Moo = 1.1, a = 10°, a combined CT = 3.0 for the four AM nozzles, and TR = 0.04 for the two south-firing ACM 
nozzles. Figure 23(a) shows a denser cluster of cells about the upper tower for the predefmed-cone grid. Details of 
the grid about the ACM plumes is seen in Figs. 23(b-d). Both the baseline and predefmed-cone grids predict 
asymmetric ACM plumes. However, the asymmetry changes sides as clearly seen in Fig. 23(d) - Station A on the 
upper tower for the finer predefmed-cone grid. The asymmetry switch is apparent in comparing the adapted grid 
views on the § = 22.5° and -22.5° planes (see Figs. 23(b-c)). The impact of the asymmetry present between the two 
ACMs is also evident in the Mach contour data at the Station B cross-section (see Figs. 23(e-f)) which occurs 
downstream of the AMs. 

Figure 24 presents the coefficient data for the OVERFLOW and USM3D solutions at M* = 1.1, -10 < a < 10°, a 
combined CT = 3.0 for the four AM nozzles, and TR = 0.04 for the two south-firing ACM nozzles. USMD data 
from both the baseline and predefmed-cone base grids of Fig. 22 are presented. As seen in Fig. 24(a), USM3D 
underpredicts C A when compared to OVERFLOW. C N and C m data have a better comparison at a < 0°. The USM3D 
C A , C N , and C m data for both grids show little variation between the baseline grid and the finer predefmed-cone grid 
(excepting the a = -5° case). 

The lateral-directional coefficients of Fig. 24(b) have been normalized by dividing each coefficient with the 
absolute value of that coefficient on the USM3D baseline grid at a = 0° which was non-zero. Using the absolute 
value at a = 0°, the asymmetry, as reflected by the sign of the coefficient, is retained. For both USM3D and 
OVERFLOW, the normalized C Y (known as A Y ) data show asymmetric flow (i.e. A Y is not equal to zero) at various 
a/CT/TR conditions. For example, at a = 10°, USM3D predicts asymmetric flow, i.e. Ay > 0.0 for the baseline grid 
which corresponds to the asymmetric jet plumes seen in the Mach contour data of Figs. 23(b-d). Also note that 
USM3D predicts a change in asymmetry between the USM3D baseline grid and the USM3D predefmed-cone grid 
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as noted by a change in the sign of A Y between the two solutions. The Mach contour data of Fig. 23(d) shows a 
corresponding switch in the jet plume asymmetry for the a = 10° case. 

The A Y data of Fig. 24(b) also show a change in angle of attack can trigger a switch in asymmetry. For example, 
OVERFLOW and USM3D both predict a change in asymmetry between the a = 0° and a = 5° conditions. Fig. 25 
presents the Mach contour data at the Station A cross-section on the upper tower for both USM3D grids at all five 
angles of attack. The USM3D predefmed-cone grid solution predicts a switch in the large jet plume from the left to 
the right side at a change in angle of attack of 0° to 5° that corresponds to a change in sign for A Y . However, the jet 
plume does not switch sides between a = 0° and a = 5° on the baseline grid. In fact at a = 0° the baseline and 
predefmed-cone grid solutions have nearly equal values in A Y while having the large jet plume appearing on 
opposites sides. Figure 26 presents a bar graph of A Y for each component of the nose-to-base geometry from both 
grids’ solutions at a = 0°. A Y for the total nose-to-base geometry is also shown on the far right in Fig. 26. The nose 
and upper tower components have equal but opposite values of A Y between the baseline and predefmed-cone grids. 
However, that trend is not necessarily the case for the components downstream of the upper tower. 

Figure 27 presents the coefficient convergence history for C A , C N , and C m for the baseline grid and the finer 
predefmed-cone grid at a = 0°. The baseline grid has a large oscillation in the coefficient history but the average of 
the coefficient over several oscillations is level and thus the solution would be considered converged. By using the 
finer predefmed-cone grid on the upper tower and thus a finer resolution for the ACM plumes, this coefficient 
oscillation is significantly reduced. Both grids give similar predictions for the averaged coefficients. 

USM3D predicted asymmetric flow for combinations of powered ACM/ AM conditions on a near symmetric 
configuration supporting asymmetric predictions by OVERFLOW. These results were used to adjust the 
aerodynamic uncertainties for powered LAV configurations. 

VIII. Concluding Remarks 

The NASA Constellation program has relied heavily on Computational Fluid Dynamics (CFD) simulations for 
generating aerodynamic databases and design loads. The NASA Tetrahedral Unstructured Software System 
(TetrUSS) tool set has served in a supporting role within the Constellation program by investigating various 
aerodynamic sensitivities and supplementing the aero database for the Orion CEV Aerosciences Project (CAP). 

The fast-paced environment of the program has driven the development of largely automated process for 
generating USM3D solutions. The flow solution process is split into 3 stages: the USM3D input stage, the grid 
generation stage, and the convergence stage. The USM3D input strategy evolved from several guideline studies and 
extensive applications. The definition of jet boundary conditions for the Abort Motor and Attitude Control Motor 
nozzles were derived by an isentropic extrapolation of nozzle exit conditions, provided by CAP, to a plane with the 
nozzle slight downstream of its throat. Grid generation evolved from relatively simple uniform grids to complex 
grids with cell clustering driven by flow structures in the USM3D solutions. The solution adaptive gridding was 
accomplished by combining the use of the TetrUSS GridTool code with the CraftTech CRISP CFD® grid adaption 
code. Solution convergence was largely based on the history of the aerodynamic coefficients. Quantitative 
characteristics of the evolving behavior for these coefficients as the flow solution develops provided a numerical 
approach to determining convergence. The adaptive gridding and solution convergence evaluation were automated 
by using scripts. Post-processing of the USM3D solution was also automated and provided the user various plots to 
aid in determining the solution’s viability. 

The TetrUSS system has provided over 2000 solutions on multiple variations of the Crew Module (CM) and the 
Launch Abort Vehicle (LAV) configurations in thirteen different tasks. Four sample applications illustrate the range 
of calculations provided by the TetrUSS system on these types of geometries. Guideline studies, such as that on the 
Apollo CM, provided direction on USM3D input and suitability of various grid approaches. By computing large 
data sets, such as that on the ID AT CM, USM3D provided corroborating data to be used in defining the aero 
database. Through comparisons with wind tunnel data such as those conducted using the 19-AA LAV wind tunnel 
model, USM3D provided not only sting interference corrections but also insight to the effects of sting interference. 
USM3D computations utilizing the automated grid generation process, such as those conducted on the ID AT CM 
configuration and the powered 606F LAV configuration, provide the ability to investigate sensitive flow features by 
variations in grid. The Orion aero-database used the USM3D data to build the aerodynamic uncertainties for both 
the CM and LAV configurations. 
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Appendix - Method for Prescribing Jet Boundary Conditions 


The methodology for defining the LAV jet boundary conditions in ideal-gas flow solvers presented below was 
provided via Excel spreadsheet by Mr. James Greathouse at NASA Johnson Spaceflight Center. The jet conditions 
are provided for by wind tunnel flow conditions or flight LAV trajectory data. A general assumption of one- 
dimensional isentropic flow is applied to determine the flow conditions at the nozzle exit plane, and then transfer 
them to the BC plane located just downstream of the jet throat. Figure A1 shows a schematic with the exit plane and 
BC plane identified. 

When using flight trajectory data, the plenum chamber has a non-air gas and thus non-air gas properties. Thus a y 
correction needs to be applied after the exit-plane conditions have been determined for the non-air jet gas. After the 
y correction, an “imaginary throat for the air nozzle” is determined using the corrected exit plane flow and the 
known exit plane area ratio. The exit plane conditions corrected for air can now be transferred to the BC plane. 

The complete jet BC calculations with a y correction are given here. For completeness, the relations for one- 
dimensional isentropic flow through a nozzle are given here: 



(Al) 

(A2) 

(A3) 

(A4) 


1. Prescribe initial conditions and parameters. 

a. Freestream conditions: 


Tqq , Poo ? Qoo ? S re f , and aoo (y air Poo / Poo) 


y air = 1.4, MW air = 28.996, R air = R universa i / MW air = 1716.47 

b. Jet plenum gas properties (for hot jet calculations): 

For hot-jet conditions: y jet =1.25, MW jet = 26.922, R jet = Runiversai / MW jet = 1846.78, T to t,jet 

For wind tunnel conditions: yj et = 1.4, MW jet = 28.996, R jet = Runiversai / MW^ = 1716.47, T to tjet 

c. Geometric areas for nozzle exit, throat, and BC plane: 

A e , A t , and A BC 

d. Assume A = A t and compute choked nozzle area ratio for the non-air jet flow: 

A e / A* = A e / A t 

2. Expand plenum chamber conditions to nozzle exit plane using jet gas properties. 

a. Using y jet , adjust nozzle exit Mach number, M e ,j e t , in Eqn. Al to match the choked area ratio with the 
geometric area ratio: 

(A e / A )j e t A e / At 

=> This defines M e je t . 

b. Substitute y je t and M e je t into Eqns. A2-A4 to compute isentropic relations at nozzle exit: 

(Ptot / Pe)jet (T tot / T e )jet (P tot / Pe)jet 

3. Calculate nozzle exit plane conditions assuming jet gas properties. 
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a. Define a target nozzle thrust coefficient: 

CT target = (target thrust) / (qoo S ref ) 

b. Adjust plenum chamber pressure, p to t in thrust equation: 

CT = [Ptot/(Ptot/Pe) jet * A e(i + 7jet M e,jet) * cos (25 <, )]/(q 00 S lef ) 

until CT matches CT ta rget • The “cos(25 )” term sets the CT in the body axial direction. 
=> This determines p tot . 

c. Calculate the total density, p toh using the gas equation: 

Ptot Ptot ! (Rjet T to t ) 

note: T to t defined in Step 1(b) 

d. Complete the nozzle exit conditions from known quantities; 


Pejet = Ptot / (p tot / pe)jet 

T e ,jet = T to t / (T to t / T e )jet 
Pejet — Pejet / (Rjet T e jet) 

V ejet — M e j et ( Yj et Rjet T e jet) 

4. Switch from jet gas to air and recalculate the chamber conditions. The general assumption is to maintain the 
thrust. Since thrust is a function of exit-plane pressure and jet momentum, the values of nozzle exit pressure, 
density, and velocity will remain unchanged. The primary impact of switching gas from jet gas to air is a change 
in the nozzle exit temperature. This directly impacts the speed of sound and exit Mach number, which in turn 
alters the values of the isentropic relations in Eqns. A2-A4 and the total conditions. 

a. Adjust nozzle exit jet temperature and Mach number for air: 

Te,air — Pejet / (Rair Pejet ) 

Me,air = V eJe t / (YairRair T e>air ) 1/2 

b. Assign remaining jet properties to nozzle exit as air: 


Pe,air Pejet 


C. 


Pe,air Pejet 

Ve,air — Vejet 

Perform sanity check of thrust by comparing new CT air with CT targ et: 


CT,: 


= [Pe,air * A e( J + y a ir M e,air) * COS ( 25 °)]/(q»S ref ) 

Compute new isentropic relations at nozzle exit using Eqns. A2-A4 with y a ir and M e;air : 

(Ptot / Pe )air (Ttot / Te )air (ptot / Pe )air 

Compute for total jet conditions with assumed air: 


Ptot, air Pe,air (Ptot / Pe )air 
Ttot, air = T e a i r (T to t / T e ) a i r 
Ptot, air — Pe,air / (ptot / Pe )air 

5. Calculate an “imaginary throat for the air nozzle”, and then A/A* at the BC plane just down stream of the 
throat. 

a. Use Eqn. A1 to compute new (A e /A*) air using y a ir and M e5 ai r . 

b. Compute new choked throat area A . 
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c. Calculate A BC /A* air 

A air - A e / (A e / A ) a i r 

d. Calculate the new choked nozzle area ratio, A B c/A a ir- 

6. Propagate the exit plane conditions to the BC plane. 

a. Using y air , adjust Mach number on the BC plane, M BC , a i r , in Eqn. A1 to match the new choked area ratio with 
the geometric area ratio (A BC / A air ). 

=> This defines M BC , a i r . 

b. Substitute y a ir and M BC , a i r into isentropic relations Eqns. A2-A4 to compute BC-plane ratios: 

(ptot / pBC )air (Ttot / TbC )air (ptot / pBC )air 

c. Compute the BC-plane properties: 

PBC = Ptot, air / ( ptot / pBC )air 
T B C = T to t, air / (Ttot / T B c )air 

pBC = P tot, air( Ptot / pBC )air 

V B C = MBc(YairRairT B c) 1/2 

d. Compute USM3D non-dimensional BC values using freestream conditions: 

PUSM3D = pBC / Poo 
VuSM3D = V B c / aoo 
PUSM3D = pBC / (Yair poo) 

7. Apply jet boundary conditions. 

a. Directly apply Pusm 3 d and Pusm 3 d across BC plane 

b. Smoothly fan the velocity vectors from the center of the BC plane until parallel to the nozzle walls. Maintain 
nozzle axis component of V US m 3 d constant across plane. 
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Table 1. List of USM3D studies. 


Task ID 

Geometry 

Title/Objective 

Guideline 1 

Apollo CM 

CEV/CM CFD Simulation Guidelines for USM3D Navier-Stokes Solver 

Develop simulation guidelines for computations on the CEV CM 

AR-10 

CEV CM 

USM3D Solutions for CEV CM Database 

Generate data for the CEV CM aerodynamic database 

Guideline 2 

Apollo LAS 

Unpowered LAS Simulation Guidelines for USM3D Navier-Stokes Solver 

Develop simulation guidelines for computations on the unpowered LAV 

AR-26 

604-Mod6 LAV 

Assessment of 604-Mod6 LAV w/o Canards or Plume Effect Using USM3D 

Generate data for an assessment database of the 604-Mod6 LAV 

AR-41 

19-AA LAV 

USM3D Analysis of the 19-AA LAV WTT 

Generate data to assess the sting influence on the 19-AA wind tunnel data 

AR-60 

605-054 LAV 

USM3D Reynolds Number Sensitivity Analysis Using the 19-AA LAV 

Generate data to assess Reynolds number effects on the LAV aerodynamics 

AR-76 

606-068 LAV 

ACM USM3D Analysis for the 606-068 LAV in Both 59-AA and Free- 
flight Configurations 

Generate data for the 59-AA LAV wind tunnel model geometry at both wind 
tunnel and flight Reynolds number conditions with the ACMs firing 

Guideline 3 

ALAS1 lr3b 

LAV with AM Jets On Simulation Guidelines for USM3D Navier-Stokes 
Solver 

Develop simulation guidelines for computations on LAV with AMs firing 

AR-100 

ALAS1 lr3b 

USM3D Flight-to-WTT Reynolds Number Scaling for 60- AA ALAS 
Including Abort Motor Jets 

Generate data for the 60-AA LAV wind tunnel model at both wind tunnel 
and flight Reynolds number conditions with the AMs firing 

AR-105 

606F LAV 

USM3D Analysis of AM Thrust Offset 

Generate data to evaluate the jet effects of the AM plume on the 606F LAV 
as a function of thrust offset 

AR-112 

606F LAV 

USM3D Analysis of AM+ACM Jet Interaction 

Generate data to evaluate the jet effects of the combination of AM and ACM 
plumes on the 606F LAV 

AR-122 

ID AT CM 

USM3D ID AT CM Cases for the v0.6 Aero Database 

Generate data for the ID AT CM aerodynamic database version 0.60 

AR-128 

26-AA LAV 

USM3D Analysis of the 26-AA WTT 

Generate data for the 26-AA LAV wind tunnel model at conditions where 
OVERFLOW did not match wind tunnel test data 
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Table 2. Information for the USM3D tasks. 


Task ID 

Geometry 

Moo 

a, deg. 

P 

deg. 

Ren 

xlO 6 

TR 

CT 

Ncell 
x 10 6 

Date 

Predefined Wake Grids 

Guideline 1 

Apollo CM 

0.7- 1.2 

0-180 

0 

8 

- 

- 

6 

3/06 

AR-10 

CEV CM 

0.7 -2.5 

0-180 

0 

10- 18 

- 

- 

6 

8/06 

Base Grids 

Guideline 2 

Apollo LAS 

0.5- 1.2 

-20-0 

0 

13-48 

- 

0 

20 

9/06 

AR-26 

604-Mod6 LAV 





- 

0 

20 

11/06 

AR-41 

605-054 LAV 

0.5 

0-180 

0 

3 

- 

0 

14 

4/07 

AR-60 

605-054 LAV 

0.5 

0, 10 

0 

3,48 

- 

0 

14 

6/07 

Solution Adapted Grids 

AR-76 

606-068 LAV 

0.3 - 0.7 

-5-20 

0, 10 

29-67 

0-0.14 

0 

27 

11/07 

Guideline 3 

ALASllr3b 

0.7- 1.2 

0- -8 


2.5 

- 

0,4 

29 

1/09 

AR-100 

ALASllr3b 

0.7 -2.5 

0- -8 

0 

2.5 - 80 

- 

0- 16 

29 

5/09 

AR-105 

606F LAV 

0.5 -3.0 

-10-10 

0 

47-83 

- 

0-4 

45 

8/09 

AR-112 

606F LAV 

0.9- 1.6 

-10-10 

0 

47-83 

0-0.1 

0-4 

58 

6/10 

AR-122 

IDAT CM 

0.5 - 8 

0-180 

0 

2.9 - 22 

- 

- 

25 

10/10 

AR-128 

26-AA LAV 

0.5- 1.3 

-15-15 

0 

2.9 

- 

1.5-6 

43 

4/11 


Table 3. Steps of the flow solution process. 


Step 

Action 

USM3D Input Stage 

1 

Create input files for USM3D: 

a. boundary condition input. 

b. jet inflow boundary condition input. 

c. flow and geometry input. 

Grid Generation Stage 

2 

Generate the baseline grid using VGRID40. 

3 

Optimize baseline grid using the CraftTech CRISP CFD® code. 

4 

Run USM3D on baseline grid for 3,000 iterations in 1 st order and then 10,000 
iterations in 2 nd order. 

5 

Perform grid adaptation using CraftTech CRISP CFD® code: 

a. adapt to gradients of p, u, v, w, and p. 

b. constrain minimum spacing to 1 inch (full scale). 

6 

Interpolate baseline solution onto adapted grid for a solution restart in USM3D. 

Convergence Stage 

7 

Restart USM3D on adapted grid and run in 1st order mode for 1,500 iterations. 

8 

Restart USM3D in 2 nd order mode and do the following: 

a. simultaneously monitor USM3D output files to check convergence criteria at 
regular timed intervals. 

b. stop solution when the first 3 convergence criteria are met to examine solution. 

9 

User examines the solution and decides from the following: 

a. extend the solution by repeating Step 8. 

b. accept the solution and begin post-processing. 
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Table 4. Convergence Criteria. 


Criterion 

Definition 

i 

Variation of C L , C D , C m , C N , and C A less than 1% over last 1000 iterations 

2 

Average of last 1000 iterations within 1% of average of preceding 1000 iterations 

3 

Maximum delta of coefficients < 0.004 over last 1000 iterations 

4 

L2-norm solution residual of all flow variables drop at least 2.5-orders of magnitude 

5 

Maximum delta of the two iteration intervals averages < 0.002 
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Service Module 


Crew Module (CM) 

\ 

Launch Abort Tower (LAT) 


Spacecraft Adapter 


Figure 1. Major components of the Orion CEV. 



(a) Side view of the Apollo CM configuration with the apex cover. 



(b) Oblique view of the ID AT CM. (c) Close-up on the nose of the ID AT CM. 

Figure 2. Views of two CM configurations. 


19 

American Institute of Aeronautics and Astronautics 



(a) Top view of the 605-054 LAV configuration. 




+X Facing 

Shoulder 

Groove 


(c) Close-up on the shoulder groove of the 606F LAV configuration. 
Figure 3. Views of two LAV configurations. 
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(d) Bottom view of the upper tower with the ACMs of the 606F LAV configuration. 



(e) Close-up view on the AMs of the 606F LAV configuration. 
Figure 3. Concluded. 



Figure 4. Sketch of predefined wake grid generation using h-refinement method. 
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Step 1 : Base Grid 


Step 2: Refine in Sphere and Cone 1 Regions 



Step 3: Refine in Cone 2 Region 


Step 4: Refine in Cone 3 Region 




(a) Grids at each step of predefined wake grid generation for a = 0°. 



(b) Predefined wake grid for a = 40°. 



(c) Predefined wake grid for a = 150°. 

Figure 5. Grids for the Apollo CM configuration at different angles of attack. 
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Base Grid 


Mach Contour Data Only 



I 


(a) Base grid and Mach contour data for the preliminary solution. 



Adapted Grid 


Mach Contour Data and Adapted Grid 




(b) Adapted grid (with and without the Mach contour data) based on the preliminary solution. 


Figure 6. Grid and Mach contour data (on the (j) = 45° plane) from the preliminary solution on the base 
grid for the 606F LAV configuration at M^ = 1.1, a = 5°, CT = 4.0, and TR = 0.04. 



Figure 7. Adapted grid (on the cj) = 45° plane) for the 606F LAV configuration at M^ = 1.1, a = 5°, CT = 4, 
and TR = 0.04. 
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Grid Adaption at I = 3 ,000 



Grid Adaption at I = 15,000 
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Figure 8. C A convergence history for two solutions obtained with different adapted grids on the 606F LAV 
configuration at = 1.1, a = 0°, CT = 3.0, and TR = 0.04. 


l st -Order Solution at I = 3,000 
Mach Contour Data on the Base Grid 


2 nd -Order Solution at I = 15,000 
Mach Contour Data on the Base Grid 




Adapted Grid Based on Solution at I = 3,000 Adapted Grid Based on Solution at I = 15,000 



(a) Close-up view on the <() = 45° plane. 

Figure 9. Mach contour data and adapted grids from the l st -order and 2 nd -order solutions for the 606F 
LAV configuration at = 1.1, a = 0°, CT = 3.0, and TR = 0.04. 
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l st -Order Solution at I = 3,000 
Mach Contour Data on the Base Grid 


2 nd -Order Solution at I =15,000 
Mach Contour Data on the Base Grid 




Adapted Grid Based on Solution at I = 3,000 


Adapted Grid Based on Solution at I = 15,000 




(b) Close-up view on the (|) = 22.5° plane. 


l st -Order Solution at I = 3,000 2 nd -Order Solution at I = 15,000 

Mach Contour Data on the Base Grid Mach Contour Data on the Base Grid 



Adapted Grid Based on Solution at I = 3,000 Adapted Grid Based on Solution at I = 15,000 



(c) Close-up view on the Station A plane. 
Figure 9. Concluded. 
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: 



Iterations 

(a) Convergence history 



Figure 10. Solution summary chart for the unpowered ALASllr3b configuration at = 1.3, a = -8°, and 
Re D = 2.5 x 10 6 at I = 9,398. 



Figure 11. Solution summary chart for the unpowered ALASllr3b configuration at M* = 1.3, a = -8°, 
and Re D = 2.5 x 10 6 at I = 11,219. 
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0 2000 4000 6000 8000 10000 

Iterations 



0 2000 4000 6000 8000 10000 

Iterations 


(c) Small oscillation convergence. 


(d) Irregular oscillation convergence 



0 10000 20000 30000 40000 50000 60000 


Iterations 

(e) Curve break convergence. 

Figure 12. Types of coefficient history convergence. 
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USM3D and Wind Tunnel Data Comparison 
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Figure 13. Comparison of USM3D and wind tunnel aerodynamic coefficient data for the Apollo CM 
configuration at = 0.9 and 0° < a <_180°. 
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Baseline Base Grid - 6,159,833 Cells 


Adapted Baseline Grid - 20,449,989 Cells 




(a) Baseline base and adapted grids. 


Finer Base Grid - 9,163,946 Cells 



Adapted Finer Grid - 32,521,627 Cells 



(b) Finer base and adapted grids. 

Figure 15. Base and adapted grids on the (|) = 0° plane for the ID AT CM configuration at = 0.5 and 
a = 150°. 


ID AT CM Geometry - Full Coefficients 
M =0.5 



Database 

— — □- — 

OVERFLOW 

o 

USM3D: Baseline Grid 

▼ 

USM3D: Baseline Grid, Time Accurate 

O 

USM3D: Finer Grid 

▲ 

USM3D: Finer Grid, Time Accurate 




145 150 

a, deg. 

Figure 16. Comparison of USM3D, OVERFLOW, and database aerodynamic coefficient data for the 
ID AT CM configuration at M^ = 0.5 and 135° < a <_160°. 
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Mach Contour Data on (|) = 0° Plane 


Baseline Grid 


Finer Grid 




Figure 17. Mach and C P contour data for the steady-state solutions on the ID AT CM configuration at 
Moo = 0.5 and a = 150°. 
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Figure 18. Coefficient convergence histories for time-accurate solutions on the ID AT CM configuration at 
Moo = 0.5 and a = 150°. 
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(a) Mach contour data on the (j) = 0° plane. 



(b) Front view of the surface C P contour data. 

Figure 19. Mach and surface C P contour data from the time-accurate solution on the finer grid for the 
ID AT CM configuration at M^ = 0.5 and a = 150°. 
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Free Air Configuration - 0° < a < 1 80° 



LAV with 15° Sting - 0° < a < 30° 



LAV with 50° Sting - 30° < a < 70° 



LAV with 162.5° Sting- 150°<a<180° 



LAV with 77.5° Sting - 60° < a < 95° 




Figure 20 Grids and a-ranges for the different model/sting configurations of the 19-AA LAV wind 
tunnel model. 
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Figure 21. Comparison of USM3D, wind tunnel, and 
model at = 0.9 and 0° < a < 180°. 


WT Data: Run 24, 15° Sting 
WT Data: Run 35, 15° Sting 
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• USM3D: Free Air 

■ USM3D: 15° Sting 

O USM3D: 50° Sting 

> USM3D: 77.5° Sting 

< USM3D: 112.5° Sting 

A USM3D: 162.5° Sting 

-* 604 Mod6 Database 



(-database data for the 19-AA LAV wind tunnel 


Baseline Grid 

Full View on 



Close-up View on 



Figure 22. Cross-section views of the baseline and 
configuration. 


Predefined-Cone Grid 

d) = 45° Plane 



the (j) = 22.5° Plane 



predefined-cone base grids for the 606F LAV 
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Baseline Grid 


Mach Contour Data 


Predefined-Cone Grid 



(a) Full view on the (j) = 45° plane. 



Adapted Grids 



(b) Close-up view on the (|) = 22.5° plane. 

Figure 23. Grids and Mach contour data at various cross-sections for the 606F LAV configuration at 
Moo = 1.1, a = 10°, CT = 3.0, and TR = 0.04. 
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Baseline Grid 


Mach Contour Data 


Predefined-Cone Grid 



Adapted Grids 



(c) Close-up view on the <|) = -22.5° plane. 




(d) View on the Station A plane. 
Figure 23. Continued. 
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Baseline Grid 


Mach Contour Data 


Predefined-Cone Grid 



Adapted Grids 



(e) Close-up view on the (j) = 45° plane. 


Baseline Grid . . . ~ A „ Predefined-Cone Grid 

Mach Contour Data 


Adapted 


(f) View on the Station B plane. 
Figure 23. Concluded. 
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606F LAV Configuration - Full Coefficients 
= 1.1, a Varies 

ACM TR = 0.04, AM CT = 3.0 


— # — 

USM3D, Baseline Grid 

— ■ — 

USM3D, Predefined-cone Grid 

♦ 

OVERFLOW 




(a) C A , C N , and C m data. 



Figure 24. Comparison of USM3D and OVERFLOW data for the 606F LAV configuration at = 1.1, 
-10° < a < 10°, CT = 3.0, and TR = 0.04. 
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Baseline Grid 


Predefined-cone Grid 


a = -10 



Figure 25. USM3D Mach contour data for two solutions on different grids at Station A on the 606F LAV 
configuration at M* = 1.1, -10° < a < 10°, CT = 3.0, and TR = 0.04. 
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Figure 26. USM3D CA data for the total configuration and each component of the 606F LAV config- 
uration at Moo = 1.1, a = 0.0°, CT = 3.0, and TR = 0.04. 
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(a) Baseline grid. (b) Predefined-cone grid. 

Figure 27. Convergence histories for the baseline and predefined-cone grids for the 606F LAV config- 
uration at Moo = 1.1, a = 0.0°, CT = 3.0, and TR = 0.04. 
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Figure Al. Schematic of nozzle jet boundary condition prescription. 
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